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Abstract. I argue that measurements of Au+Au collisions at 20, 130 and 200 GeV of 
the centrality dependence of the mean pt together with pt and net-charge fluctuations 
reflect the approach to local thermal equilibrium. 
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>• ' 1. Introduction 

OO ' 

O . Fluctuations of the net transverse momentum reported by the CERES, NA49, PHENIX 



and STAR experiments exhibit substantial dynamic contributions pQ. In particular, 



O . 

preliminary STAR data for Au+Au collisions at 20, 130 and 200 GeV energies show 
that p t fluctuations increase as centrality increases l2J. However, data from PHENIX 
and STAR exhibit a similar increase in the mean transverse momentum (p t ), a quantity 
unaffected by fluctuations [3 II]- In I ask whether the approach to local thermal 
equilibrium can explain the similar centrality dependence of (pt) and pt fluctuations. 
Here I address new data from STAR and PHENIX, including a new and, perhaps, 



related effect in the centrality dependence of net-charge fluctuations (2). 

Dynamic fluctuations are generally determined from the measured fluctuations by 
subtracting the statistical value expected, e.g., in equilibrium 0. For particles of 
momenta pi and P2, dynamic multiplicity fluctuations are characterized by 

Raa = ^ =jNjtJ d ^ d P 2 r ^ ^ C 1 ) 

where (• • •) is the event average. This quantity depends only on the two-body correlation 
function r(pi,p 2 ) = AT(p 1 ,p 2 ) — AT(p 1 ) A/(p 2 ). It is obtained from the multiplicity 
variance by subtracting its Poisson value (N), and divided by (N) 2 to minimized the 
effect of experimental efficiency jB]. For dynamic p t fluctuations one similarly finds 

(SptiSpti) = J dp l dp 2 j^^jr ^tife, (2) 

where 5pu = pu — {p*}; STAR measures this observable. The observable measured by 
PHENIX satisfies F pt pa N (5ptiSpt2) /2cr 2 when dynamic fluctuations are small compared 
to statistical fluctuations a 2 = (p 2 ) — (pt) 2 - 
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2. Thermalization 

Thermalization occurs as scattering drives the phase space distribution within a small 
fluid cell toward a Boltzmann distribution that varies in space through the temperature 
T(x, t). The time scale for this process is the relaxation time z/ _1 . In contrast, density 
differences between fluid cells must be dispersed by transport from cell to cell. The 
time needed for diffusion to disperse a dense fluid "clump" of size L ~ (iVnl/n)^ 1 is 
td ~ uL 2 /v^ h , where t> t h ~ 1 is the thermal speed of particles. This time can be much 
larger than v~ x for a sufficiently large clumps. Global equilibrium, in which the system is 
uniform, can be only obtained for t 3> td- However, the rapid expansion of the collision 
system prevents inhomogeneity from being dispersed prior to freeze out. 

Dynamic fluctuations depend on the number of independent particle "sources." 
This number changes as the system evolves. Initially, these sources are the independent 
strings formed as the nuclei collide. The system is highly correlated along the string, 
but is initially uncorrelated in the transverse plane. As local equilibration proceeds, the 
clumps become the sources. The fluctuations at this stage depend on number of clumps 
as determined by the clump size, i.e., the correlation length in the fluid. Dynamic 
fluctuations would eventually vanish if the system reaches global equilibrium, where 
statistical fluctuations are determined by the total number of particles. 

3. Mean p t and its Fluctuations 

In jH] I use the Boltzmann transport equation to show that thermalization alters the 
average transverse momentum following 

{pt) = {Pt)oS+{ Pt )e(l-S), (3) 

where S = e _A/ " is the probability that a particle escapes the collision volume without 
scattering. The initial value (pt) is determined by the particle production mechanism. 
If the number of collisions M is small, S ~ 1 — jV implies the random- walk-like increase 
of (p t ) relative to (pt) - For a longer-lived system, energy conservation limits (p t ) to a 
local equilibrium value (pt) e fixed by the temperature. 

As centrality is increased, the system lifetime increases, eventually to a point where 
local equilibrium is reached. Correspondingly, the survival probability 5* in decreases 
with increasing centrality. The average pt peaks for impact parameters near the point 
where equilibrium is established. The behavior in events at centralities beyond that 
point depends on how the subsequent hydrodynamic evolution changes {pt) e as the 
lifetime increases. Systems formed in the most central collisions can experience a cooling 
that reduces (|HJ) with proper time r as (pt) e oc r~ 7 [5]. 

Thermalization can explain the behavior in fig. ^ The survival probability is 
S = exp{— z/(r)dr} « (to/tf)", where v = (av ve \)n{r) for a the scattering 
cross section, v re \ the relative velocity, and to,f the formation and freeze out times. 
Longitudinal expansion implies n(r) oc r -1 , yielding the power law with a = z/ r o- 
To fit the measured centrality dependence, I assume a = 4 and 7 = 0.15 in central 
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Figure 1. Mean p t (left) and dynamic p t fluctuations (right) with data 0]. 
The dashed curve in the right figure is a fit to the PHENIX data alone, while the solid 
curve is a fit to the STAR data. Centrality is determined by the number of participants 
relative to the b=0 value. 



collisions, and parameterize S(N part ) by taking a oc N art and rp — r oc iVp ar t, where 
^part is the number of participants. I take the same a for all species, as appropriate for 
parton scattering. 

Dynamic fluctuations depend on two-body correlations and, correspondingly, are 
quadratic in the survival probability. In [3] I add Langevin noise terms to the Boltzmann 
equation to describe the fluctuations of the phase space distribution. A simple limit is 
obtained when the initial correlations are independent of those near local equilibrium, 
(SptiSp t 2) = (SptiSp t 2) S 2 + {5pti5pt2) e (l — S) 2 ; this form was used in Alternatively, 
if the initial correlations are not far from the local equilibrium value, I find 

(6 P nS Pt2 ) = (5 Ptl 5p t2 ) S 2 + (SpaSpaUl - S 2 ). (4) 

Here I use (|4*]) : which provides somewhat better agreement with the latest STAR data 
in the peripheral region where thermalization is incomplete and my model assumptions 
most applicable. As before, the initial quantity (5pt\5p t2 )o is determined by the particle 
production mechanism, while (Spti5pt2) e describes the system near local equilibrium. 

To estimate (5pu5p t2 ) for nuclear collisions, I apply the wounded nucleon model 
to describe the soft production that dominates (pt) and (SptiSpt2), to find 

2(5pti5p t2 )p P fl + R pp \ 

{SpnSpa} ° = iv, [muJ (5) 

The pre-factor is expected because (J2J) measures relative fluctuations and, therefore, 
should scale as N~^ t . The term in parentheses accounts for the normalization of (j2J) to 
{N(N— 1)); Raa scales as N~^ rt [6j. ISR measurements imply (SptiSpt2)p P /(pt)l P ~ 0.015. 
I use HIJING to estimate R pp and Raa, thus building in resonance and jet fluctuations. 

Near local equilibrium, spatial correlations occur because the fluid is inhomogeneous 
- it is more likely to find particles near a dense clump. These spatial correlations 
fully determine the momentum correlations since the distribution at each point is 
thermal. The mean p t at each point is proportional to the temperature T(x), so 
that (SptiSp t 2)e ~ / r ( x i, x 2)^( x i)^^( x 2)j where r(x 1 ,x 2 ) is the spatial correlation 
function, 5T = T — T, and T is a density-weighted average. I take n (oc T 3 ) and r to 
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Figure 2. Dynamic pt (left) and net-charge fluctuations (right) with STAR data 2 . 
The curves and data are for three different beam energies, 20 GeV, 130 GeV and 200 
GeV. Centrality is determined by the number of participants relative to the b=0 value. 



be Gaussian with the transverse widths R t and £ t , respectively the system radius and 
correlation length. In [1 I obtain 

(SpnS P n) e = F { f^ (6) 

J- + tiAA 

where Raa is given by (0) [5 . The dimensionless quantity F depends on £t/Rt- I 

1/2 

compute F assuming R t oc N' t ] F = 0.046 for £t/Rt = 1/6. I again use HIJING to 
estimate Raa- 

Calculations in figs. ^ and EJ illustrate the common effect of thermalization on one- 
body and two-body pt observables. The solid curves in all figures are fit to STAR 
fluctuation data and (p t ) data, while the dashed curve shows a fit to the PHENIX data 
alone. New to this work are the STAR p t and net charge fluctuation data in fig. |2l To 
compute the energy dependence of fluctuations, I take a oc N and Raa oc N" 1 for the 
measured multiplicity N. Deviations in the most central collisions at the highest beam 
energies may result from radial flow or, perhaps, jets pQ. Nevertheless, I stress that jets 
- and quenching effects - are incorporated in my calculations via the HIJING Raa m 
(0). HIJING without rescattering does not describe this data. 

Net charge fluctuations in fig. [21 are characterized by v dyn = R ++ + R — 2R + ~ 
for R ab given by ((TJ) [H]. I therefore expect v dyn to satisfy v dyn = v Q S 2 + u e (l — S 2 ). I 
use STAR pp data to fix uq and, following [Hj, compute v e /v assuming a longitudinal 
narrowing of the correlation length seen in balance function data pQ . The agreement of 
in-progress calculations with preliminary data [2] is encouraging. 
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